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ABSTRACT 

Context. In classical T Tauri stars. X-rays are produced by two plasma components: a hot low-density plasma, with frequent flaring 
activity, and a high-density lower temperature plasma. The former is coronal plasma related to the stellar magnetic activity. The latter 
component, never observed in non-accreting stars, could be plasma heated by the shock formed by the accretion process. However its 
nature is still being debated. 

Aims. Our aim is to probe the soft X-ray emission from the high-density plasma component in classical T Tauri stars to check whether 
this is plasma heated in the accretion shock or whether it is coronal plasma. 

Methods. High-resolution X-ray spectroscopy allows us to measure individual line fluxes. We analyze X-ray spectra of the classical 
T Tauri stars MP Muscae and TW Hydrae. Our aim is to evaluate line ratios to search for optical depth effects, which are expected 
in the accretion-driven scenario. We also derive the plasma emission measure distributions EMD, to investigate whether and how the 
EMD of accreting and non accreting young stars differ. The results are compared to those obtained for the non-accreting weak-line 
T Tauri star- TWA 5. 

Results. We find evidence of resonance scattering in the strongest lines of MP Mus, supporting the idea that soft X-rays are produced 
by plasma heated in the accretion shock. We also find that the EMD of MP Mus has two peaks: a cool peak at temperatures expected 
for plasma heated in the accretion shock, and a hot peak typical of coronal plasma. The shape of the EMD of MP Mus appears to 
be the superposition of the EMD of a pure coronal source, like TWA 5, and an EMD alike that of TW Hydrae, which is instead 
dominated by shock-heated plasma. 

Key words, stars: atmospheres - stars: coronae - stars: individual: MP Muscae, TW Hydrae, TWA 5 - stars: pre-main sequence - 
techniques: spectroscopic - X-rays: stars 

1. Introduction ature Tmax = 3yumHVQ/(16fe), where vq is the infall velocity. 

^, . , _ „ . , Because of the high pre-shock velocity, the infalling material 

Classical T Tauri stars (CTTS) are young low-mass stars, still ^^^^^^^ temperatures of a few MK, and hence it emits X-rays, 

surrounded by a circumste lar disk from which they accrete ma- ^^p^^^j ^^j^^^ ^^^^ ^^^^.^^.^^^ ^^^^ ^.^^^ .^^.^^^^ ^^^^ 



terial. According to a widely accepted model, they have a strong ^j^^ accretion-dri ven X-ray luminosity should be comparable to 
H magnetic field that regulates the accretion process disrupting the ^^^^^^^ (lGullbrind ri994.). Considering typical inferred 

circumstellai- disk, loading material of the inner part of the disk, stream cross-section al area (< 5% of the stellar surface, e.g. 

and guid ing it in a free fall along i ts flux tubes toward Aecgn- |Calvet & Gullbring|[T998h . velocity (~ 300 - 5 00 km s'), and 

talstar ( lUchida & Shibata| [IM i Bertout et al. |[198i |Koenigl| '^^ss accretion rate (~ lO--^ - 10'^ vr'^ e.g. lCxullbrinJet al.l 

' f ■ , . . ,. ,. Il998l) . it can be inferred that the plasma heated in the accretion 

A characteristic feature of young stai's is strong X-ray emis- ^j^^^,^ ^j^^^j^ j^^^^ ^^^^.^.^^ > lO'^ cm-^, i.e. at least one or- 

sion that IS traditionally ascribed to magnetic activity in their magnitude higher than coronal plasma density. Hence, in 

coronae. Similar to their niore evolved siblings the diskless ^^^^.^^.^^ 

weak-hne T Tauri stars (WTTS), CTTS display high X-ray lu- j^. temperatures of a few MK. 

minosities and frequent flaring activity. The typical temperature _ . . . ^ ^rr^r^o • • 

of their coronal plasma is ~ 10 -20MK,orev en higher during , ^o suinmarize. X-ray emission from CTTS can originate 

sti-ong flai-es (e.g. 50- 100 MK, Getman et al.ll2008h . ^""^ two different plasma components, plasma heated m the ac- 

c »u 1 • » f ■ »u cretion shock and coronal plasma. The former, because of its 

From a theoretical point of view, the accretion process can , , , , , • , ,V , , , 

also produce significant X-ray emission on CTTS. Material ac- l"^^'" temperatures, should dominate the softer X- ray band (e.g 

f ,u ■ . 11 J- 1 u 1 f nnn E < IkeV in the case of the CTTS TW Hya, Gunther et al. 

creting from the circumstellar disk reaches velocities of ~ 300 - ottt^ ■ ■ ^ ^ — \ . . , 

500kms-i. A shock forms at the base of the accretion col- ™- ^hile the harder X-ray emission, E > 1 keV, should be 

umn because of the impact with the stellai- atmosphere. This P^duced almost entirely by coronal plasma. 

shock heats up flie accreting material to a maximum temper- Recently, high-resolution X-ray spectra of a few CTTS 

enabled measurement of individual emission lines sensitive 

Sen^/ Oj^nw? re^Mfijfj fo: C. Argirofli, e-mail: argi@astropa.unipa.it to plasma density (i.e. He-like tiiplets), and hence seai'ches 
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for evidence of accretion-driven X-rays. The density of the 
plasma at T ~ 2 - 4MK can be inferred from the Ovii 
and Neix triplet lines (at E » 0.6 and 0.9 keV, respec- 
tively). All but one of the CTTS for which the Ovii triplet 
lines were detected showed cool plasma with high density, 
»e > 10" cm "^ iKastner et al. 2002; Stelzer & Schmitt 2004; 
Schmitt etal.' '2005; Guntheretal. 2006; Argi roffi et all |2007|; 
Gudel et al. 2007; Robrade & Schmitt 2007). In contrast, the 
cool quiescent plasma of active stellar coron ae is always dom- 
inated by low de nsities («e ~ 10'° cm'-^. iNess et al.l l2002t 
iTesta et al.ll2004al) . This basic difference suggests that the high- 
density cool plasma in CTTS is not coronal plasma but plasma 
heated in accretion shocks. One complication to this argument is 
that mass accretion rates derived from assuming a very high ef- 
ficiency of conversion of accretion energy into X-rays tend to be 
an order of magn itude or so lower than rates deri v ed using other 
metho ds (e.g. iDr^ 120051; ISchmitt et al] 120051; iGiinther et all 
I2007h . 

The idea of accretion-driven X-rays from CTTS is su- 
perficially supported by a soft X-ray excess found in high- 
resolution X-ray spectra of CTTS with respect to similar spec- 



tra of WTTS by Telleschi et al. (2007D and 
(120071) . However, Gudel et al.. {,200% and 



Giidel & Telleschi 
Giidel & Telleschii 



(I2007h noted that this soft X-ray excess is significantly lower 
than that predicted by simple models of X-ray emission from ac- 
cretion shocks. Moreover, the soft excess scales with total stel- 
lar X-ray luminosity, and hence is related at least partially in 
so me way with the stel l ar ma gnetic activity. iGiidel etal] (l2007h 
and lGtidel & Telleschii (|2007|) suggested that the CTTS soft X- 
rays could be produced by infalling material loaded into coronal 
structures. 

The properties of the X-ray emitting plasma in CTTS and 
in WTTS have also been investigated using CCD X-ray spec- 
tra of large stellar samples. These studies, however, commonly 
covered the 0.5 - 8.0 keV energy band in which the coro- 
nal component dominates. The main results are that CTTS are 
on average less lumin ous in the X-ray band than WTTS (e.g. 
iFlaccomio et" aL"2003'). and that X-ray emitting plasma o f CTTS 
is on the avera ge hotter than that of WTTS (Neuhaus er et al.l 
[l99llPreibisch et al.i2005h . CTTS and WTTS therefore do have 
difl^erent coronal characteristics, suggesting that the accretion 
process can affect coronal properties to some extent. 

Numerical simulations have confirmed that the accretion 
process can produce significant X-rays: Giinther et al. (2007) de- 
rived stationary 1-D models of the shock in an accretion column; 
ISacco et al] (l2008a l) improved those results by performing 1-D 
hydrodynamical (HD) simulations of the accretion shock, in- 
cluding the stellar atmosphere and taking in to account time vari- 
ability. Assuming optically thin emission, ISacco et al.l (l2008al) 
showed that, even for low accretion rates, the amount of X- 
rays produced in the accretion shock is comparable to the typ- 
ical X-ray luminosity of CTTS (Lx ~ lO^^ergs^' for M ~ 
lO-i^Moyr-i). 

Several aspects of the nature of the high-density cool plasma 
component observed in CTTS are still debated. In particular, 
definitive evidence that it is material heated in the accretion 
shock is still lacking^Moreover, while the simple "photospheric 
burial" model of llJrakei (i2005b suggests that under some circum- 
stances a large fraction of the shock X-rays can be absorbed and 
reprocessed by the photosphere, there are currently no detailed 
quantitative models explaining why the X-ray luminosities, pre- 
dicted on the basis of 1-D HD simulation results and mass ac- 
cretion rates inferred from observations at other wavelengths, 
are universally much higher than observed. 



Understanding the link between accretion and X-rays would 
also allow more accurate characterization of the coronal com- 
ponent of the X-ray emission from CTTS. This could help 
in understanding how accretion changes coronal activity, and 
which other parameters determine the coronal activity level 
in PMS stars, whose X-ray luminosity cannot simply be ex- 
plaine d in terms of a Rossby dynamo number dPreibisch et al.l 
20051) as is largely th e case for active main sequence stars (e.g. 
Pizzolato et al.ll2003l) . 

To address the above issues, we performed a detailed study 
of the high-resolution X-ray spectra of two nearby CTTS: 
TW Hya and MP Mus. In particular we investigated: 

- optical depth effects in their soft X-ray emission; 

- the emission measure distribution (EMD) of the X-ray emit- 
ting plasma. 

Optical depth effects probe the nature of the high-density 
cool plasma component: we show that, if the emitting plasma 
is located in the accretion shock, some emission lines should 
have non-negligible optical depth; in contrast these Unes should 
be optically thin if the plasma is located in coronal structures. 
We also investigate how the EMD can help in recognizing coro- 
nal and accretion plasma components, which should have dif- 
ferent average temperatures. We compare the EMD of accret- 
ing and non-accreting young stars, and compare these observed 
EMD with that pred icted on the basis of the HD shock model of 
(ISacco et al.ll2008ah . 



2. Targets 

We selected the two CTTS MP Mus and TW Hya for 
this study. These stars suffer moderate interstellar absorption 
(A^H < 10^' cm"^), and good S /N spectra, gathered with XMM- 
Newton/RGS (high spectral resolution and large effective area 
in the soft X-ray band) are available. While high-resolution X- 
ray spectra of TW Hya have also been obtained by Chandra, we 
chose to analyze the XMM-Newton/RGS spectra so as to have 
more uniform data for the two stars, enabling more ready com- 
parison of the derived results. 

TW Hya is a 0.7 M© CTTS, located at 56 pc. From its 
membership in the eponymous TW Hya Association ( TWA) , 
its age is estimated t o be ~ lOMyr ( Kastner et al.l 119971) . 
iMuzerolle et al.l (l2000h estimated a mass accretion rate of ~ 
5 X 10~"'Mcr)yr~ ' based on its Hg profi l e and UV flux. 
lAlencar & Batalhal (12002') and 'Batal ha et all (|2002|) derived a 
higher value, 1 - 5 x 10 '^MQyr based on UV data, in 
agreement with the mass accretion rate in ferred from the Ha 
line width by iJavawardhana et aTl (l2006h . The X-ray emis- 
sion of TW Hya shows clear evidence of high-density plasma 



Table 1. Log of flie XMM-NewtonfRGS observations of MP Mus 
and TW Hya. 



Instrument Start 


Exposure 




(UT) 


(ks) 


MP Mus (ObsId 0406030101) 


RGSl 


2006 Aug 19 09:30:48 


98.7 


RGS2 


2006 Aug 19 09:30:56 


98.8 


TW Hya (ObsId 01 12880201 ) 


RGSl 


2001 Jul 9 05:51:39 


29.0 


RGS2 


2001 Jul 9 05:51:39 


28.0 
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Fig. 1. RGS spectra of MP Mus and TW Hya. Each spectrum was obtained by adding RGS 1 and RGS2 spectra, rebinned by a factor 
2. Labels identify some of the strongest lines. 

Table 2. Strongest RGS lines of MP Mus and TW Hya. 



^pred 


Ion 


log r* 

max 




MP Mus 
Flux'' 


EMiy 


^obs 


TWHya 
Flux'' 


EMD'' 


8.42 


Mg XII Mg XII 


7.00 








8.39 


4.1 ± 4.1 




10.24 


Ne X Ne X 


6.80 








10.28 


7.6 ± 5.2 


* 


11.00 


Neix 


6.60 








11.01 


10.8 ± 7.2 


* 


11.54 


Neix 


6.60 








11.55 


23.6 ± 16.7 


* 


12.13 


Ne X Ne X Fe xvii 


6.80 


12.14 


22.2 ± 3.4 


* 


12.15 


58.6 ± 9.4 


* 


12.28 


Fe xxi Fe xvii 


7.00 


12.30 


5.5 ± 2.6 


* 








12.85 


Fe XX Fe xx Fe xx Fe xx 


7.00 


12.86 


3.4 ± 2.3 










13.45 


Ne IX Fe xix 


6.60 


13.46 


21.1 ± 3.7 




13.46 


244.7 + 18.5 


* 


13.55 


Ne IX Fe xix 


6.60 


13.55 


8.6 ± 3.1 




13.56 


162.4 + 16.4 




13.70 


Neix 


6.60 


13.73 


8.8 ± 2.8 




13.71 


91.7 + 29.8 




14.21 


Fe XVIII Fe xviii 


6.90 


14.19 


3.3 ± 1.7 


* 








14.82 


O VIII viii 


6.50 








14.84 


10.5 + 4.8 


* 


15.01 


Fe XVII 


6.70 


15.02 


11.8 ± 2.9 


* 


15.05 


26.5 ± 5.7 


* 


15.18 


O VIII VIII Fe XIX Fe xvii 


6.50 


15.21 


7.6 ± 2.5 


+ 


15.21 


23.6 ± 7.7 


* 


16.01 


O VIII VIII Fe XVIII Fe xviii 


6.50 


16.02 


22.6 ± 6.3 


* 


16.03 


43.2 ± 15.5 


* 


16.78 


Fe XVII 


6.70 


16.78 


8.8 ± 1.8 




16.79 


13.3 ± 6.8 


* 


17.05 


Fe XVII Fe xvii 


6.70 


17.08 


17.2 ± 3.0 


* 


17.06 


24.3 + 7.6 


* 


17.77 


Ovii 


6.40 








17.80 


5.7 + 6.6 


* 


18.63 


Ovii 


6.30 


18.65 


8.7 ± 2.6 




18.65 


21.6 + 5.8 


* 


18.97 


O VIII vm 


6.50 


18.98 


65.6 ± 4.7 




18.98 


292.7 ± 13.3 


* 


19.83 


N VII N VII 


6.40 








19.85 


4.3 ± 4.4 


* 


20.91 


N VII N VII 


6.30 








20.95 


14.3 ± 9.8 


* 


21.60 


Ovii 


6.30 


21.60 


30.2 ± 4.0 


* 


21.62 


191.7 ± 17.5 


* 


21.80 


Ovii 


6.30 


21.81 


27.9 ± 7.9 




21.83 


104.6 ± 14.3 




22.10 


Ovii 


6.30 


22.10 


8.0 ± 2.8 




22.12 


13.8 ± 8.7 




24.78 


N VII N vn 


6.30 


24.80 


10.5 ± 2.8 


* 


24.79 


118.2 ± 13.4 


* 


28.47 


C VI C VI 


6.20 


28.49 


3.8 ± 1.3 


* 








28.79 


Nvi 


6.20 








28.82 


25.2 + 5.2 


* 


29.08 


Nvi 


6.10 








29.10 


11.5 ± 5.0 




33.73 


C VI C VI 


6.10 


33.74 


16.4 ± 5.2 


* 


33.75 


86.5 ± 13.8 


* 



" Predicted (APED database) and observed wavelengths (A). * Temperature (K) of maximum emissivity. 
with uncertainties at the 68% confidence level. Lines selected for the EMD reconstruction. 



Observed line fluxes (10 * ph s 



(we ~ 10^-^ cm ^) at l ow temperatures jKastner et al.l 120021; 
IStelzer & Schmitai2004l) . 

MP Mus is a CTTS with spectral type KlIVe. lMamaiek et alj 
(120021) identified it as a member of the Lower Centaurus Crux 
(LCC) Association, and determined a distance of 86 pc using the 



moving cluster method and the LCC convergent point. They in- 
ferred for MP Mus a mass of 1.1 - 1.2 Mq from photometry 
and isoch rone fitting to different evolutionary tracks. Recently 
[Torres et a l. (2008) suggested that MP Mus is likely a member 
of the younger e Cha Association. The authors, using their con- 
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vergent m ethod iteratively w ith their list of candidates of the as- 
sociation (iTorres et al.l2006h . determined a distance to the star of 
103 pc. Note that both these distance estimates, 86 and 103 pc, 
depend on correct identification of the parent stellar associa- 
tion. In this work we assume a distance for MP Mus of 86 pc. 
Adopting the longer distance would lead to higher X-ray fluxes 
by only 30% or so and would be of no consequence for the pur- 
pose of this study. Like TW Hya, MP Mus also ex hibits evidence 
of coo l plasma at high density («e » 5 x 1 " cm ^ , lArgiroffi et al.l 
l2007h . While it i s certain that MP Mus is a CT TS (Ho- equivalent 
width is ~ 47 A, iGregorio-Hetem et al. Ill992h . Its mass accretion 
rate is not known. The accretion rate is fundamental for studying 
the relation between accretion and soft X-ray emission, and we 
evaluate it in Sect.|4]using an optical spectrum of MP Mus taken 
with FEROS at La Silla Observatory. 

We compare the results obtained for TW Hya and MP Mus 
with those obtained for TWA 5, a member of the TWA with 
similar age to TW Hya. TWA 5 is a quadruple system with 
three similar components of 0.5 Mq and a brown dwarf. Ho- 
equivalent width (~ 1 A. G regorio-Hetem et al.l 1 99 2) indicates 
that TWA 5 is a WTTS, even if some margin al accretion is 
still p resent at least in one of the components jMohantv et al.l 
We selected TWA 5 for the comparison because: 1) 
its EMD was derived from an XMM-Newton/RGS observation 
( I ArgirofR et al .1120051) with the same method we adopted in this 
work for MP Mus and TW H ya; 2) low plasma de nsity was in- 
ferred from the Ovii triplet (lArgiroffi et al.l l2005h . suggesting 
that its X-ray emission is entirely produced by the stellar corona. 



Table 3. EMD and abundances of MP Mus and TW Hya. 





MP Mus 


TW Hya 


logr(K) 


log EMD (cm-^) 


log EMD (cm-^) 


6.00 
6.20 
6.40 
6.60 
6.80 
7.00 
7.20 


5i.76!;;{^ 

52.24!JJ:42 

5 3. 03 

53 00+° °^ 

jj.uu Q20 

52 76+""' 

-"^■'"-0.08 

52 89+" ''' 

53 13+"" 


51 79+"-25 

— U.Uo 

51 ■74^0^ 
52.64!g:?4 

CO 00+0.05 

J^-^^-o.oe 
51 41+"i2 

-'^•^^-0.09 

51 63+""^ 
ji.ojm3 

51 69+""^ 


Elem. 


Ax/AxQ 


Ax/Axe 


C 

N 

O 
Ne 
Mg 

Fe 


0-28^!^ 

18+009 

46+014 
= Fe 

09+""" 


70+" '^'^ 

"•'"-0.17 
o oT+0.45 
^•^■-'-0.44 

96+"-'' 
7 66+1-3" 

' -""-0.94 

2 30+'-^4 

43+0.00 



Abundances are in solar units jAsplund et al.ll2005b . EMD val- 
ues have been derived assuming that MP Mus and TW Hya are 
located at 86 and 56 pc, respectively. 



3. RGS data analysis 

TW Hya was observed on July 2001 for 30 ks, and MP Mus 
on August 2006 for llOks. Detai l s of the d ata processing 
can be found in Stelzer & Schmitt' ('2004') and lArgiroffi et al.l 
( l2007h . respectively. In Table [1] we report the log of the XMM- 
Newton/RGS observations of the two stars. For MP Mus we dis- 
carded time segments affected by high background count rates, 
obtaining a net exposure of ~ lOOks over the selected good time 
intervals. 

Background-subtracted spectra of MP Mus and TW Hya 
contain ~ 5800 and ~ 6500 net counts, respectively. To increase 
the S/N ratio of the data we added the RGSl and RGS2 spectra 
and rebinned by a factor 2. The resulting RGS spectra are shown 
in Fig.[T] 

We based the spectral analysis (line identification, line 
emissivity functions G(T), line oscillator strengths /) on 
the Astr ophysical Plasma Emission Database (APED VI. 3, 
ISmith et al. 2001), as implemented in the PlNTofALE V2.0 
software package (Kashyap & Drake 2000). Spectral line fluxes 
were measured within PlNTofALE by fitting the observed lines 
with a Lorentzian profiles. The results are listed in Table|2] 

3. 1 . EMD reconstruction 

Having defined a temperature grid, T,, with binsize AT, the EMD 
of the X-ray emitting plasma is computed as the amount of 
emission measure, rig rin AV, of plasma with temperature rang- 
ing between T,- - Ar/2 and T, + AT/2. We derived the EMD 
from the line flux measure ments using the Markov-C hain Monte 
Carlo (MCMC) method of Kash vap & Drakdd 19981) . We consid- 
ered all the measured lines whose intensity depends only on the 
plasma temperature, and not on the plasma density: hence, we 
discarded the intercombination and forbidden lines of the O vii 



and Ne ix He-like triplets. The lines used for the EMD recon- 
struction are marked in Table |2] 

We reconstructed the plasma EMD over a logarithmic tem- 
perature grid ranging between log T(K) = 6.0 and 7.2, with a bin 
size of 0.2. This choice was guided by the set of formation tem- 
peratures of the measured lines (see log r^ax in Table |2]i. This 
range includes the expected temperature of plasma heated in an 
accretion shock. On the other hand, coronal plasma may also 
have temperatures higher than log T = 7.2, so the derived EMD 
does not take very hot coronal components into account. 

The derived EMD for MP Mus and TW Hya are reported 
in Table[3] Simultaneously with the EMD reconstruction, we de- 
rived also the relative abundances of elements represented by the 
selected lines. The absolute abundances, also reported in Table[3] 
were fixed for both stars by matching the predicted and observed 
spectra. 

3.2. Optical deptli effects 

Assuming that the X-ray emitting plasma has a geometrical 
depth / (in cm) along the line of sight, free electrons with den- 
sity «e and temperature T ^e (in K), the n the optical depth of an 
emission line is given bv (lActonlll978l) : 

r=1.16xl0-i/ter'^^^«e/ (1) 
\ Te y nz nil "e 

where A is the line wavelength (in A), / the line oscillator 
strength, niz is the atomic mass of the element (in amu), «z,, 
the density of element Z with ionization level /, nz the density of 
element Z, ne the hydrogen density (all the densities are in units 
of cm"-') 

For increasing density n^, and/or increasing source dimen- 
sion I, T increases. A non-negligible optical depth {t ~ 1), be- 
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cause of resonance scattering, first occurs in lines with large os- 
cillator strengths / and ground state lower levels. Photons pro- 
duced in these transitions are absorbed again by ions of the same 
species in the ground state. These photons are then re-emitted in 
different, random directions. When the optical depth is of the or- 
der of unity, resonance scattering only affects the emergent spec- 
trum if the source does not have a spherical symmetry. 

In the scenario of accretion-driven X-ray emission, the 
material heated in the accretion shock is entirely located in 
a small compact volume at the footpoint of the accretion 
stream. The high-density cool plasma component observed in 
CTTS has density of 10" - 10'^ cm"^ and linear dimension of 
/ ~ 10^ - 10'" cm (inferred from measured values of «e and 
£M.lKastner et al.ll2002l:[Stelzer & Schm ittll2004tlAr giroffi et all 
I2007r iRobrade & Schmitlll2007l) . Considering these values for 
«e and /, the optical depth of the strongest emission lines pro- 
duced by this shock-heated plasma (T ~ 1-3 MK) should be 
non-negligible. As an example, the line center optical depth, t, 
of the Ovii resonance line at 2 1.60 A is t ~ 10 for a den- 
sity Hf. - I X 10" cm"-', rizj/nz - 0.5 (corresponding to a 
temperature T = 2MK), mz/hh = 8.5 x lO""*, rin/n^ = 0.83, 
mz = 16.0amu, and / = 10^ cm. 

On the other hand, if this high-density plasma component 
observed in CTTS is contained in many separate coronal struc- 
tures, the average length / covered by photons before escaping 
from the plasma should be significantly smaller, and hence the 
optical depth negligible. In fact, extensive investigation of coro- 
nal spectra for optical depth effects has shown coronae to be ef- 
fectively optically-thin in almost all cases (e.g. Ness et al. 2003; 
iTesta et aiTl2004bl l2007h . with only a few notable exceptions 
(see below). 

Optical depth effects in a given line can be investigated by 
considering, as a reference, another line produced by the same 
element in t he same ionization stage, but with a smaller oscilla- 
tor strength (iNess et al.ll2003l: iTesta et al.ll2004bl: iMatranga et"aLl 



Table 4. Line ratios sensitive to optical depth. 



120051: iTesta et al.l 20071) . In the optically-thin case, the ratio be- 
tween two such lines is dictated purely by atomic parameters and 
has only a weak dependence on the plasma temperature. Instead, 
when the gas is not optic ally- thin, the line with larger oscillator 
strength can be quenched relative to the weaker line, causing a 
discrepancy between the observed and predicted optically-thin 
ratios. 

In the high-resolution X-ray spectra gathered with Chandra 
and XMM-Newton, the strongest emission lines with the largest 
oscillator strengths are the Lye lines of O vm and Ne x, the reso- 
nance lines of Ovii and Neix, and t he Fe xvii lin e at 15.01 A. 
Consideririg Cha ndmlWEYGS, data, iTesta et al.l [2004b ) and 
iTesta et akl (l2007h found significant resonance scattering in the 
Lyo'/Lyyg line rati os of O viii and N e x o f the ac tive stars II Peg 
and IM Peg. MatrangaetaD ( |2005|) and iRose et al. (2008), an- 
alyzing RGS data of AB Dor and EV Lac, found that the ra- 
tio between the Fexvii lines at 15.01 and 16.78 A changes be- 
tween the flaring and quiescent phases, indicating opacity in the 
15.01 Aline. 

We searched for optical depth effects in the X-ray emission 
of MP Mus and TW Hya, using TWA 5 as a non-accreting com- 
parison. We investigated three line ratios: 

- the ratio of the Ovii lines at 21. 60 A and 18.63 A, pro- 
duced by transitions to the ground level (n = 1) from the 
shells n - 2 and n = 3 respectively; the resonance line at 
21.60 A has a large oscillator strength (/ = 0.69), while the 
line at 18.63 A has a small oscillator strength (/ = 0.15); 



Star 


A^H (cm-2) 


Ovil(r) 


Ovin(LvQ-) 


Fexvn(15.0) 








MP Mus 


5 X 10^° 


4.0+1.5 


3.2 + 1.0 


1.3 + 0.4 


TW Hya 


2 X 10^0 


9.4 + 2.8 


7.0 + 2.6 


1.9+1.1 


TWA 5 


3 X 10^° 


8+4 


11+4 


1.4 + 0.4 


Procyon 


1.1 X 10"* 


7.9 + 1.1 




1.2 + 0.4 


AT Mic 


1 X 10"* 


6.2 ± 2.6 




1.8 + 0.6 


5CrB 








2.4 + 0.4 


YZCMi 


1.2 X 10"* 






1.2 + 0.4 


44 Boo 


1 X 10'^ 






2.39 + 0.28 


a Cen 




9.7+1.0 






Capella 


1.8 X 10'^ 






1.64 + 0.12 


AB Dor 


1.7 X lO'"^ 


8.9+1.6 




1.85+0.10 



,Active stellar 



Ratios are derived from line : 
measurements are taken from: Raassen et a!. (2002) for Procyon; 
Raassenetal. (2003) for AT Mic; Gondoin (2005) for 6 CrB; 
Raassen pt aT.I i?mih for YZ CMi '; GondoinI ilOQ^Jo^M Boo; 



Audard et alJ (120011) for Capella; iSanz-Forcada et all (l2003i) for 
AB Dor. 



- the ratio between the O vin Lyor at 18.97 (a doublet with os- 
cillator strengths / = 0.54 and / = 0.27) and the O vm Ly^S 
at 16.01 A (a doublet with / = 0.10 and / = 0.05); 

- the ratio between the Fe xvii lines at 15.01 A (large oscillator 
strength, / - 2.73) and 16.78 A (small oscillator strength, 
/ = 0.11). 

This line set allowed us to probe the optical depth of the plasma 
at temperatures ranging from 2 to 5 MK. We did not consider the 
Ne X LycK and Ly/3 line ratio since, for MP Mus, we detected only 
the Lycc line. 

The observed ratios of MP Mus, TW Hya, and TWA 5 are 
listed in Table |4] and plotted in Fig. |2] These ratios were ob- 
tained correcting the observed fluxes (Table |2]i for the interstel- 
lar absorption, assuming the A^h values derived from the analysis 
of EPIC spectra (Argirofli et al. 2007; Stelzer & Schmitt 200^ 
I Argiroffi et al.l 120051) . The correction factors due to interstellar 
absorption range between 0.94 and 1.16. We associated each 
flux ratio of each star to an average temperature. That average 
temperature was obtained assuming, as a weighting function, the 
product of the EMD and emissivity functions, G{T), of the line in 
the numerator. The predicted flux ratios as a function of plasma 
isothermal temperature are also illustrated in Fig.|2] 

All the lines involved in the analyzed ratios are free from 
significant blending, except for the Ovm Ly/3 line, located at 
16.01 A. This line is blended with two Fexvm lines (16.00 and 
16. 07 A), whose co ntributions are small but non-negligible (see, 
e.g. iTesta et aTl 1200 7, for a detailed discussion and treatment): 
the EMD and abundances derived for the X-ray emitting plasma 
of MP Mus and TW Hya (see Sect. 13. It indicate that the Fe xviii 
contribution to the observed fluxes is about ~ 20% and ~ 6%, re- 
spectively. To take this into account, we included these contribu- 
tions in the calculation of the predicted flux ratios. The predicted 
ratio of the O vm Lya and Ly/3 lines is then weakly dependent 
on the model abundances and hence differs sUghtly for different 
stars. 
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Fig. 2. Line photon flux ratios sensitive to optical depth. Filled 
symbols mark the observed line ratios, corrected for interstellar 
absorption, for the three stars MP Mus, TW Hya, and TWA 5. 
Each ratio is associated with the average temperature obtained 
assuming the function EMD(T) x G(T), as weighting function. 
The temperature associated to the observed ratios for the active 
stars (open symbols) selected is the temperature corresponding 
to the maximum value of the emissivity function G{T) of the 
21.60 and 15.01 A lines. Predicted values are computed assum- 
ing optically thin emission from an isothermal plasma. Upper 
panel: photon flux ratio of the O vii lines at 21.60 and 18.63 A. 
Medium panel: photon flux ratio of the O viii Lya vs. Ly^S lines. 
Due to a blend with Fe lines, predicted values differ for the three 
stars because of the different relative abundances. The solid line 
indicates the ratio predicted for MP Mus, the dotted line (barely 
visible below the solid line) for TW Hya, and the dashed line 
for TWA 5. Lower panel: photon flux ratio of the Fe xvii lines at 
15.01 and 16.78 A. 



For comparison, we considered also other coronal sources 
for which line fluxes and A^h values, obtained from XMM- 
Newton/RGS data, were pubhshed. These stars and their O vii 
and Fe xvii line flux ratios, corrected for the interstellar absorp- 
tion, are listed in Tableland are included in Fig.|2] We did not 
consider the O viii line ratio since its analysis requires knowl- 
edge of model abundances for each star 

MP Mus shows in all cases an observed ratio significantly 
lower than the predicted one, over the entire range of plau- 
sible temperature. The discrepancy between the observed and 
predicted ratios for MP Mus are 2.8cr, 4.7cr, and 2.3cr for the 
O VII, O VIII, and Fe xvii ratios, respectively. Conversely, the ob- 
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Fig. 3. PN (black) and MOSl (red) spectra of M P Mus with su- 
perimp osed the 3-T model (solid line) derived bv lArgiroffi et al.l 
( l2007l) indicating low absorption (A^h = 5x10^° cm"^), and the 
3-T model (dashed line) obtained by fitting the EPIC spectra 
with an absorbing column fixed at 5 x 10^' cm"^. 

served ratios for TW Hya are perfectly compatible with those 
predicted for optically-thin emission. The O vii and O viii ratios 
for all the other stars considered are compatible with the pre- 
dicted optically-thin values. Conversely some stars, other than 
MP Mus, show a Fe xvii line ratio discrepant with the optically 
thin case. The interpretation of the Fe xvii ratio is more contro- 
versial: while coronal plasma might have non-negligible optical 
depth in the 15.01 A line (i .e . the Sun and AB Dor, Sab a et al.l 
119991: iMatranga et aTll2005l) . iNess et alj (l2003l) found that sev- 
eral stars, in a large sample, exhibit ratios somewhat discrepant 
from the predicted optically-thin value, as in our case, but they 
concluded that the interpretation of this result in terms of opacity 
effects is not convincing. 

Summarizing, we detected a significant intensity deficit in 
lines with the largest oscillator strength in the spectrum of 
MP Mus: the O vii resonance line at 2 1 .60 A, the O viii Lyo' line 
at 18.98 A, and the Fe xvii line at 15.01 A. These lines are mainly 
formed by the plasma at a temperature of 2 - 8 MK. We conclude 
that this emission is not optically-thin and that resonance scat- 
tering has quenched the strongest lines. 

We explored whether the anomalous oxygen line ratios ob- 
served in MP Mus could be explained by an absorbing column 
A^H higher than that assumed (that however could not explain 
the Fexvii line ratio). To explain the observed Ovii and Oviii 
line ratios for MP Mus, a hydrogen column ^ 5 x 10^' cm"^ is 
required. This value is higher, by a factor 10, than the value con- 
straine d from the analysis of the EPIC spectra bv lArgiroffi et al.l 
(1200 7'). and adopted for the line ratio analysis. In Fig.[3]we show 
the PN and MOSl spectra of MP Mus (MOS2 is omitted for 
clarity). We also plotted the 3-r m odel with the inferre d ab- 
sorption A^h = (5 ± 1.8) X 102"cm-2 dArgiroffi et al.ll2007l) . and 
the 3-r model obtained by fitting the EPIC spectra with an 
absorbing column fixed at 5 x 10^' cm"^. The soft part of the 
EPIC spectra, which constrains the A^h value, is dominated by 
continuum emission. Therefore, if soft continuum and soft RGS 
lines, both contained in the 0.3 - 1.0 ke V energy range, are pro- 
duced by the same plasma component, then the observed line 
ratios of MP Mus cannot be explained in terms of photoelectric 
absorption. Instead if soft continuum and soft lines originated 
from different plasma components, affected by different absorb- 
ing columns, then the observed line ratios could be explained by 
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photoelectric absorption, rather than opacity effects in the emit- 
ting plasma. However, different absorbing columns would indi- 
cate that the highly absorbed component (the one producing the 
O VII and O viii lines) is located under the accretion stream or 
buried down in the stellar atmosphere, suggesting again that this 
plasma component originates in the accretion shock. 

4. Optical spectra and their analysis 

To gather information on the accretion status of MP Mus, we 
analyzed its optical spectrum taken from the public archive of 
La Silla Observatory (ESO). The observation was performed on 
March 18, 2005 using the echelle spectrograph FEROS at the 
2.2m telescope. The spectra cover approximately 3850-9200 A 
with a resolving power R - 48000 (AA 0.15 A, in the region 
of Ha, measured in the calibration lamp spectrum). This range 
includes all the chromospheric activity indicators, ranging from 
Ca II H & K to the calcium infrared triplet, as well as the lithium 
Une at 6708 A. The signal-to-noise ratio at 6550 Ais S /N ~ 15. 

For the reduction, we used the standard procedures in the 
IRAFQ package (bias subtraction, extraction of the scattered light 
produced by the optical system, division by a normalized flat- 
field, and wavelength calibration). After reduction, the spectrum 
of MP Mus was normalized to the continuum order by order by 
fitting a polynomial function to remove the general shape from 
the aperture spectra. We did not perform a flux calibration since 
it is not necessary for our study. 

4.1. Ha line and accretion rate 

We measured the equivalent width of the Ha line using an IDL 
procedure developed by us that enables estimation of the error 
in the measurements using both the signal-to-noise ratio and the 
spectral resolution. We obtained EW(Ha) - -41.05 + 0.09 A, 
which is similar to the value reported by Mamajek et al. (2002), 
but lower than that measured by IGregorio-Hetem et alj (I1992h . 
who found EW(Ha) = -47 A. This variation of 6 A in obser- 
vations obtained ~ lOyr apa rt is not very high if compared e.g. 
with variations observed by Sacco et al. (2008b) in cr Ori and 
A Ori. However, there is not enough time coverage to study the 
Ha variability in detail. 

We also measured the width at 10% of the peak of the Ha 
line (see Fig. |4]i, which is indica tive of the ma ss accretion rate 
for widths above 200kms"' (Natta etal] 12004 ). We obtained 
W(Ha) = 446 + 5 kms"', corresponding to a mass accretion 
rateof ^ 3 x lO^'^Moyr-i. 

5. Discussion 

5.1. Soft X-rays from siiock-lieated plasma: optical depth 
effects 

It appears that in most CTTS, significant soft X-rays are pro- 
duced by a high-density plasma component. This result has been 
interpreted as a strong indication that this plasma, at a tempera- 
ture of a few MK, is formed in the shock at the base of the ac- 
cretion stream near the stellar surface. Hence this high-density 
plasma component should be contained in a small volume at the 
base of the accretion stream, instead of being located in extended 
structures like coronal plasma. 




6560 6570 
Wavelenght (Angs) 



' IRAF is distributed by the National Optical Observatory, which is 
operated by the Association of Universities for Research in Astronomy, 
Inc., under contract with the National Science Foundation. 



Fig. 4. Ha line profile of MP Mus, illustrating the line width at 
10% of the maximum. 



We showed in Sect. 13.21 that, considering the characteristic 
volumes and densities for the high-density cool plasma compo- 
nent of CTTS, the optical depth t (eq.[Tll in some emission lines 
with large oscillator strength should be significantly larger than 
1, producing detectable opacity effects. We found clear evidence 
of opacity effects in all the examined line ratios of the soft X- 
ray spectrum of MP Mus, produced by the high-density plasma 
at temperatures of a few MK. These lines are: the resonance 
lines produced by n = 2 to n = 1 shell transitions in the O vii 
and Oviii ions (at 2 1.60 A and 18.97 A, respectively), and the 
Fexvii line at 15.01 A. All these lines are significantly weaker 
than their predicted optically -thin intensities, a result which can 
be readily explained in terms of resonance scattering. We did not 
find evidence of optical depth effects in the X-ray emission from 
TW Hya. 

We explored whether the observed opacity could be pro- 
duced by coronal plasmas. Studies on large stellar surveys 
showed that optical depth effect s from cool coronal plasma ex - 
ist but are extremely rare (.Ness et al.l 120031: iTesta et aL I l2007h . 
Quiescent plasma confined within coronal magnetic structures 
can hardly explain t he observed non-negligible optical depths 
in CTTS. In fact the iRosner et a l. (1978) model indicates loop 
semi-length of ~ lO*" - 2 x 10^ cm, for T,^^^ = 3 MK and = 
5xlO"-10'-' cm"-'. Such scale lengths are unrealistically small 
if compared to photospheric and chromospheric scale heights, 
and not large enough to account for the observed optical depth (t 
would be ~ 0.03, in the center of the O viiline at 21.60 A, assum- 
ing that the loop cross section has a radius smaller by a factor 10 
than the loop semi-length). In principle flaring plasma, instead 
of quiescent plasma, could favor opacity effects. In fact stan- 
dard models suggest that plasma density increases during flares, 
raising the opacity. O pacity effects have been searched in flaring 
coronal sources (e.g. lGiidel et al.l2004|). but thev have been de- 
tected only a couple of times dMatranga et alj|2005l : iRose et all 
2008). In those cases opacity affected Fexvn lines formed at 
higher temperatures than the O vii and O vm lines. Moreover the 
soft X-ray light curve of MP Mus ind icates that this emis sion 
does not originate from flaring plasma dArgiroffi et al.ll2007l) . 

We cannot reject altogether the hypothesis that the high- 
density cool plasma of MP Mus, showing opacity effects, origi- 
nates from coronal plasma. However the alternative explanation, 
based on the hypothesis that this plasma is located in the post 
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Fig. 5. Schematic view of the accretion stream footpoint at the 
stellar surface. The gray volume indicates the post-shock region 
containing the high-density plasma at temperatures of a few MK. 



shock region at the base of the accretion stream, naturally ex- 
plains the observed opacity effects. Therefore opacity by itself is 
not a prove of the nature of this plasma component, but strongly 
supports the shock-heated plasma hypothesis. 

We find evidence of optical depth effects for MP Mus, and 
not for TW Hya, although optical depth effects were expected 
for both stars. Moreover the higher densities of the plasma 
of TW Hya, in comparison to MP Mus, should favor opac- 
ity effects, considering that the characteristic dimension of the 
plasma volume / scales with EM^^^ n^^^^, and hence t oc In^ - 
EM^I^ n\^^ . A tentative explanation of the lack of opacity effects 
from TW Hya could be related to the fact that resonance scatter- 
ing depends on the geometry of the source and on the inclination 
angle under which the source is observed. In particular, the post- 
shock accretion hot-spot is likely to have quite different vertical 
and horizontal dimensions. 

The volume occupied by the plasma heated in the accretion 
shock is defined by the stream cross-section A and by the thick- 
ness of the hot post-shock region L (see Fig.|5]l. L is mainly de- 
termined by the post-shock velocity and b y the plasma coolin g 
time, and hence by the plasma density (e.g. lSacco et alj|2008ah . 
Estimates of both A and L can be made from the soft X-ray emis- 
sion of CTTS: observed electron densities range from 10'' cm""* 
to lO'^ cm"^, a typical value of the EM - nl x A x L respon- 
sible for the soft X-ray emission is ~ 10^^* cm"-'. Assuming an 
infall velocity of 400 km s"' and assuming that all the shock- 
heated plasma is contained in one shock area, we obtain: L ~ 
10^ - 10'' cm and A ~ 10^° - lO^^cml The difference of the 
horizontal and vertica l dimensions does not d epend on assuming 
only one shock area dRomanova et alj|2004l) . inste ad of assum- 
ing a few separated ( < 10) accretion footpoints jPonati et alj 
l2007h . 

Therefore, the post-shock region has the vertical dimension 
L significantly smaller than the other two dimensions. In Fig.|5] 
we display a schematic view of the footpoint of an accretion 
stream and the post shock region. The resonance scattering effect 
on the X-ray spectrum should be different depending on whether 
we observe the emitting plasma from the top or from the side 
(inclination angle 0i or 62 respectively). The latter is the one that 
favors optical depth effects, while in the former case these effects 
could even be negligible. The different opacity effects observed 
for MP Mus and TW Hya could be explained assuming that for 
the two stars the angle 6 of the accretion footpoints differ signif- 
icantly. TW Hya and MP Mus are viewed under quite different 
inclinations /: TW Hya is almost pole-on (/ ^ 10°, e.g. Oi et al. 
|2004 . while MP Mus should have / ^ 60°, considering that its 
optical spectrum indicates vsin/ ~ 13kms"' and that its rota- 



tional period is ~ 5 - 6 d dBatalha et aljl998l) . Assuming that the 
footpoints of the accretion streams are located at high latitudes 
on both stars (i.e. the magnetic dipole field i s almost aligned 
with the rotational axis. iRomanova et al.ll2004l) . the orientation 
of TW Hya should be similar to the top- view {0i in Fig.|5]l, while 
that of MP Mus to the side view (62)- We suggest that this im- 
portant difference in viewing angle might explain the different 
optical depth obtained from the X-ray spectra of the two accret- 
ing stars. We cannot exclude that, although the inclination angles 
i of two stars differ significantly, the two viewing angles 9 of the 
post shock region are similar (e.g. if the accretion footpoint lat- 
itude is ~ 50 - 60°). If it were the case, a different explanation 
would be required. 

5.2. Soft X-rays from shock-heated plasma: EMD 

The above results are consistent with the current prevailing view 
that accreting stars have two distinct plasma components: coro- 
nal plasma and shock-heated plasma, both contributing to differ- 
ent extents to the observed X-ray emission. To study the proper- 
ties of these two plasma components it is necessary to identify 
and disentangle them with the help of both density and tempera- 
ture diagnostics. 

The EMD can also provide information on the nature of the 
X-ray emitting plasma. In Fig.|6]we show the EMD of the two 
CTTS MP Mus and TW Hya. To understand their EMD in terms 
of the coronal and accretion-driven hot pl asma components, w e 
compared them with the EMD of TWA 5 ( Argiroffi et al.ll2005h . 
rebinned on the same log T grid used for the EMD of MP Mus 
and TW Hya. 

The EMD reconstruction is based on the assumption that 
emission is optically thin. We know that this hypothesis is not 
true in the case of MP Mus. The lines affected are the resonance 
lines of O vii and O vm which constrain the cooler part of the 
EMD. Hence, it must be noted that the true EMD of MP Mus 
at low temperatures (log T $ 6.5) is likely to be slightly higher 
than that inferred by a factor of 2 or so. 

Inspecting Fig.|6]we found three main results: 

1. MP Mus and TW Hya show a peak in their EMD at low 
temperature (log T ~ 6.4); this peak is not present in the 
EMD of TWA 5; 

2. MP Mus and TWA 5 have a strong EMD peak at high tem- 
perature (logT ~ 7.0 - 7.2), while for TW Hya that peak, 
even if marginally detected, is significantly lower compared 
to those of the other two stars; 

3. the relative strength of cool and hot peak EMD in MP Mus 
and TW Hya is significantly different. 

ISacco et akl (l2008al) performed hydrodynamic simulations, 
tuned to the MP Mus case, to study the shock formed by an ac- 
cretion stream impacting on a chromosphere and its X-ray emis- 
sion. In Fig. Q we show the EMD of MP Mus together with th e 
time-averaged EMD of the simulation by Sacco et all (l2008al) . 
The normalization of the model EMD was derived assuming that 
the soft X-ray emission is entirely produced by the shock-heated 
plasma component. This is reasonable, considering the O vii line 
triplets: any contribution from coronal p lasma, assu ming that it 
has low density, cannot exceed 20% ( Argiroffi et al.li2007b . The 
important result is that the model EMD has a pronounced peak 
at r ~ 3 MK, which is exactly the feature that we observe in 
MP Mus and TW Hya, while no cool peak is present in the EMD 
of the WTTS TWA 5. 

The cool peak that we found in the EMD of the two CTTS 
TW Hya and MP Mus explains, in terms of EMD, the soft flux 



C. Argiroffi et al.: X-ray Optical Depth in T Tauri Stars 



9 



53.5 
^ 53.0 
^ 52.5 

Q 

^ 52.0 
cn 
o 

51.5 
51.0 - 



MP Mus 
TW Hyo 
TWA 5 



:4 



6.0 



6.5 7.0 

log T (K) 



7.5 



Fig. 6. Emission measure distributions derived from the RGS 
line fluxes. 

Table 5. Mass accretion rate. 
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exce ss observed in the X-ray spectra of CTTS in th e XEST sur- 
vey (iGtidel & Telleschill2007tlTelleschi et al.l l2007'). 

These results suggest that the shape of the EMD in CTTS 
can be used to determine whether soft emission can be produced 
by accretion shocks, also in those cases in which the electron 
density provided by the O vii trip let is compa tible with coronal 
plasma (i.e. the case of T Tau, Giidel et al. 2007), or when the 
vital diagnostic intercombination and forbidden lines are of poor 
quality. 



5.3. Soft X-rays from shock-heated plasma: mass accretion 
rate 

The comparison between the EMD of MP Mus with that of the 
WTTS TWA 5, and the results obtai ned from the h ydrodynam- 
ical simulations performed by Sacco et alJ (l2008ah shows that 
the soft X-ray emission of MP Mus can be entirely explained 
in terms of the plasma heated in the accretion shock. The mass 
accretion rate provided by the simulations, 8 x 10"" Moyr"', 
agrees with th a t deriv ed on the basis of a simplified model by 
lArgiroffi et al.l (l2007l) . However the M value derived from the 
Ha line width exceeds that obtained from X-rays, Mx, by more 
than one order of magnitude. This discrepancy was noted also 
by [Drake (2005) and Gunther et al. (2007) for TW Hya. A sim- 
ilar situation arises for BP Tau, for which ma ss accretion rates 
were derived from X-ray data by .Schmitt et al.l (12005) . We list in 
Table |5] mass accretion rates derived from X-ray data for these 
three CTTS, compared to those derived from other accretion in- 
dicators (i.e. UV or Ha). In all three cases the Mx values are 
lower by a factor 10 (or even more) than the corresponding value 
derived from UV or optical data. 
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Fig. 7. EMD of MP Mus derived from the RGS line fluxes, and 
the EMD deri ved from hydrodyna mical simulations of the accre- 
tion shock bv lSacco et al. I (l2008al) . 



We note that the observed fluxes of the strongest lines were 
used to derive Mx, assuming optically thin emission. Since we 
found that the X-ray emission produced by the shock-heated 
plasma, in MP Mus, is not optically thin, the accretion rate 
of MP Mus derived from X-rays, (5 - 8 x lO^^Mgyr"', de- 
rived from the Lyg and resonance line s of Ovm and Ovii, 
lArgiroffi et all 120071; ISacco et all l2008al) . is Ukely underesti- 
mated by a factor of 2 or so — much too httle to reconcile the 
discrepancy between the M values. 

The accretion flow is likely composed of several funnels, 
each isolated from the adjacent one because the strong magnetic 
field inhibits thermal conduction and mass motion perpendicu- 
lar to its lines. Each funnel is characterized by a given density 
and infall velocity. Both the density and velocity determine the 
amount of observable X-rays produced by the accretion shock 
of each funnel. As explained below a distribution of density and 
velocity might explain the observed M discrepancy. 

Drake (2005) argued on the basis of accretion stream ram 
pressure and shock stand-off height that more dense streams 
would form smaller, less extended shocks buried too deeply in 
the photosphere to be observed from most inclination angles. X- 
rays would then be reprocessed to lower energy by the surround- 
ing photospheric gas. The observed accretion-shocked plasma 
of MP Mu s should have a stand-off height of a 10 ^ - 10^ cm 
or so (e.g. lArgiroffi et al.ll2007l: ISacco et alj|2008al) — suflicient 
to lie well above the photosphere and be visible with minimum 
absorption from the stellar atmosphere. It is possible that there 
are other streams of accretion carrying the majority of the in- 
flowing mass flux that are much more dense and that do get ob- 
scured from the line-of-sightby absorption. X-ray emission from 
TW Hya apparently does not fit in this scenario: it is produced 
by a very high- density plasma, but it suffers low absorption. To 
solve this issue. iDrakd (l200 5) suggested that the plasma density 
of TW Hya could be slightly overestimated, because of UV ra- 
diation and its influence on the He-like triplets. 

Another ef fect mig ht also contribute to the observed discrep- 
ancy. IRomanova et all (120041) and FOregorv et al.l (l2006l) infeiTed 
that accreting material arrives at the base of the accretion stream 
with a large distribution of infall velocities. Therefore assuming 
that different funnels (or different part of the same funnel) may 
have different velocities, then each funnel will produce a post 
shock region with a temperature depending only on the relevant 
velocity. In such case, funnels with high velocity would produce 
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hot post shock plasma, and, hence. X-ray emission, while fun- 
nels with reduced velocity might produce post-shock plasma of 
insufficient temperature to contribute to the X-ray emission. 

5.4. Coronal X-ray luminosity in CTTS 

The presence of coronal plasma for both TW Hya and MP Mus 
is indicated by their observed flaring activity and by the pres- 
ence of a hot plasma component detected in their EPIC spec- 
tra, and inferred by the presence of emission lines prod uced by 
hot pl asma (T ~ 10 MK, Kastner et al. 2002; Stelzer & Schmitj 
l2004t lArgiroffi et al .1 120071) . Assuming that the cool peak in the 
EMD of MP Mus and TW Hya is entirely due to accretion, 
then their coronal component can be evaluated considering only 
the hottest part of their EMD. If we take into account only 
log T(K) > 6.7, the coronal X-ray luminosity of MP Mus and 
TW Hya are 1.4 x 10^^° and 1.4 x lO^'^ergs"', respectively, in 
the 0.5 - 8.0 keV band. The coronal luminosity of TW Hya is 
significantly reduced if compared to its whole X-ray luminosity, 
8.1 X 10^'^ ergs"', while for MP Mus the corona emits ~ 80% 
of its entire X-ray luminosity. With these new estimates for the 
coronal Lx, MP Mus and TW Hya have log(Lx/iboi) = -3.4, 
and -3.8 respectively. Both of these values, but especially that 
of TW Hya, are significantly lower than previous estimates of 
their coronal to bo lometric luminosities (Ma maiek et al.l 120021; 
iKastner et al.l 11997.) . and place the two CTTS under the sat- 
uration level log(Lx /Lboi) = -3, as is usual for CTTS (e.g 
iPreibisch et al.ll2005l) . This suggests that, when the X-ray lumi- 
nosity is to be used as activity indicator, it would be preferable 
to compute it excluding the plasma components at T < 5 MK. 

6. Conclusions 

In this work we presented an analysis of the high-resolution X- 
ray spectra, obtained with XMM-Newton/RGS, of the two CTTS 
TW Hya and MP Mus. For MP Mus we detected significant res- 
onance scattering in the Oviii Lyo' and Ovn resonance lines, 
which are produced by the high-density plasma at temperatures 
of a few MK. No resonance scattering was detected in the spec- 
trum of TW Hya. This result strongly supports the hypothesis 
that this plasma is formed by the shock at the base of the accre- 
tion column which is likely viewed obliquely. The different opti- 
cal depths observed for TW Hya and MP Mus could be explained 
in terms of different viewing angle of their accretion shocks. 

We also derived the EMD for TW Hya and MP Mus, finding 
that they both show a peak at T ~ 3 - 4 MK, in addition to the 
hot peak at 10 - 20 MK typical of coronal plasma on magneti- 
cally active stars. The cool peak is perfectly described by plasma 
heated in an accretion shock ( Sacco et al..2008 a). The same peak 
is not present in the plasma of the non-accreting young star 
TWA 5. The identification of this EMD peak as due to shock- 
heated plasma allows us to assess the characteristics of the true 
coronal plasma and its luminosity. In particular, the coronal X- 
ray luminosity of TW Hya is less than 20% of its whole X-ray 
luminosity in the 0.5 - 8.0keV band. 

Soft X-ray emission can be used to compute the mass ac- 
cretion rate. We compared the mass accretion rates derived from 
X-ray data, Mx, with those obtained from UV and optical indi- 
cators, M, finding that Mx is underestimated by, at least, a factor 
10. Two possible explanations for this are that some accretion 
streams are most dense and form shocks that are located too deep 
in the stellar atmosphere to be observed, and that some plasma 
in accretion streams might not attain the full free-fall velocity 
required to form X-ray emitting shocks. 
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